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ABSTRACT 

The angular power spectrum of the Sunyaev-Zel'dovich (SZ) effect is a powerful probe of 
cosmology. It is easier to detect than individual clusters in the field, is insensitive to obser- 
vational selection effects and does not require a calibration between cluster mass and flux, 
reducing the systematic errors which dominate the cluster-counting constraints. It receives 
a dominant contribution from virialized cluster region between 20 — 40% of the virial ra- 
dius and is thus relatively insensitive to the poorly known gas physics in the cluster centre, 
1^ . such as cooling or (pre)heating. In this paper we derive a refined analytic prediction for the 

■ SZ angular power spectrum using the universal gas-density and temperature profile and the 

| dark-matter halo mass function. The predicted power spectrum has no free parameters and 

f^*) ■ fits all of the published hydrodynamic simulation results to better than a factor of two for 

2000 < I < 10000. We find that the angular power spectrum C/ scales as C; cx cr|(£lb/i) 2 
and is almost independent of all of the other cosmological parameters. This differs from the 
local cluster abundance studies, which give a relation between erg and Q m , We also com- 
pute the covariance matrix of C; using the halo model and find a good agreement relative to 
the simulations. We argue that the best constraint from the SZ power spectrum comes from 
I ~ 3000, where the sampling variance is sufficiently small and the spectrum is dominated by 
massive clusters above 10 14 h~ x Mq for which cooling, heating and details of star formation 
are not very important. We estimate how well we can determine aslfl^h) 2 ^ 7 with sampling- 
k>( \ variance-limited observations and find that for a several-square-degree survey with arcminute 

resolution one should be able to determine cr§ to within a few percent, with the remaining 
uncertainty dominated by theoretical modeling. If the recent excess of the CMB power on 
small scales reported by CBI and BIMA experiments is due to the SZ effect, then we find 
<7 8 (ft b /i/0.035) - 29 = 1-04 ± 0.12 at 95% confidence level (statistical) and with a residual 
10% systematic (theoretical) uncertainty. 

Key words: cosmology: theory — cosmic microwave background — dark matter — galaxies: 
haloes — galaxies: clusters: general — cosmological parameters 
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1 INTRODUCTION 

Abundance of dark-matter halos and its redshift evolution have 
been recognized as an important cosmological test. For example, 
the local abundance of X-ray clusters can be used to constrain 
a combination of the matter-fluctuation amplitude and the matter 
density of the universe (e.g., White, Efstathiou, & Frenk 1993; 
Viana & Liddle 1996; Eke, Cole, & Frenk 1996; Kitayama & Suto 
1997). Higher-redshift clusters can be used to break this degeneracy 
(e.g., Bahcall & Fan 1998), as well as constrain other cosmological 
parameters (Haiman, Mohr, & Holder 2001). 
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There are two major uncertainties in such analysis. One is the 
sample itself, which has to be carefully constructed to avoid any 
possible selection biases that can affect the final constraint. The 
other is the relation between observable quantities (e.g., X-ray tem- 
perature or flux) and the halo mass, which is related to the theo- 
retical models. The latter uncertainty is particularly important as 
a small error in this relation can result in a significant effect on 
the final results. For example, the X-ray flux scales as gas-density 
squared and is dominated by the central regions, where complicated 
physical processes are taking place. This means that the relation be- 
tween flux or flux-weighted temperature and the halo mass is not 
easily obtained from ab-initio numerical simulations, and one must 
resort to more empirical calibrations instead. 

The Sunyaev-Zel'dovich (SZ) effect has some advantages 
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over X-ray in this respect. The SZ effect is caused by the cosmic 
microwave background (CMB) photons scattering off hot electrons 
in the intracluster medium, and the SZ flux is proportional to the 
projected electron-gas pressure, which is proportional to a product 
of electron number density and temperature along the line of sight. 
Relative to X-ray emission this has a larger contribution from the 
outer parts of the cluster, where the missing physics such as gas 
preheating or cooling in numerical simulations is less significant. 

However, the SZ effect is very difficult to measure and cur- 
rently only detections of already-known clusters exist (Birkinshaw 
1999). It is likely that the first SZ detection in a random direction 
will not be that of a specific cluster, but of a field correlation func- 
tion, or the angular power spectrum Ci, with several low signal-to- 
noise clusters contributing to the signal. Even if individual clusters 
are too faint to be detected with high statistical significance, the 
combined effect in Ci can be statistically significant. There have 
been several observational efforts in measuring the small-scale fluc- 
tuations, and upper-limits have been obtained by ATCA (Subrah- 
manyan et al. 1998), BIMA (Holzapfel et al. 2000; Dawson et al. 
2001), SUZIE (Church et al. 1997; Ganga et al. 1997), and VLA 
(Partridge et al. 1997). Recently, CBI and BIMA experiments have 
reported statistically significant detections of the small-scale fluc- 
tuations at the level of 15 — 20 piK on arcminutes angular scales 
(Mason et al. 2002; Dawson et al. 2002), and their results may al- 
ready be providing the first detection of the SZ fluctuations (Bond 
et al. 2002). 

There are other advantages of using power spectrum in ad- 
dition to the higher statistical significance in a noisy map. One is 
that selection effects are less important than in a cluster survey, for 
which flux, surface brightness, and other selection criteria have to 
be carefully examined to avoid any biases. The other is that we do 
not have to measure the mass of individual haloes. As discussed 
above this is a dominant source of uncertainty in the local-mass- 
function determination from X-ray clusters and lack of a good cal- 
ibration can lead to a 30% error on as- Even if mass were deter- 
mined from SZ flux, for which the central regions are less impor- 
tant, one would still need to include effects such as substructure, 
projections, ellipticity etc. As we argue in this paper using the SZ 
power spectrum instead of the number-count analysis is less sen- 
sitive to selection biases and can still give statistically significant 
results. It is true that the power spectrum does not use all of the 
information available, specially if one obtains the information on 
cluster redshifts. However, the SZ power spectrum provides quite 
powerful constraints on cosmological parameters already by itself 
and the small increase in statistical errors is more than compen- 
sated by the decrease in the systematic uncertainties. Thus initial 
experimental efforts should focus on the power spectrum as a more 
reliable method to determine the cosmological parameters. 

There are of course other effects that contribute to fluctuations 
in this range of wavelengths. The SZ effect is the dominant con- 
tribution for / > 2000 — 4000, depending on its amplitude. The 
primary CMB dominates on the larger scales. Although the pri- 
mary CMB confuses the SZ effect, we can reduce or eliminate it by 
observing at 217 GHz at which the SZ effect vanishes. Many of the 
current generation CMB experiments (such as CBI) operate at low 
frequencies, where the frequency dependence of CMB and the SZ 
effect is similar, and the two cannot be distinguished very well on 
this basis. Galactic emission (synchrotron, dust, or free-free) does 
not have much power on the small angular scales, and has a dif- 
ferent frequency dependence from CMB. They are not likely to be 
a major source of contamination. We ignore the kinetic SZ effect, 



as it is at least an order of magnitude smaller in power than the 
thermal SZ effect. 

Point sources have a lot of small-scale power and are poten- 
tially more problematic contaminant. While we can easily identify 
bright (more than several-a) point sources, we cannot identify the 
fainter ones, which will contaminate the SZ effect. This is particu- 
larly problematic for narrow frequency coverage experiments such 
as CBI, as we cannot use different frequency dependence of the 
point sources to separate the two components. We can still try to 
use opposite signatures between SZ and point sources as a way to 
distinguish them; the SZ effect is negative in flux at low frequen- 
cies, while point sources are positive. This results in a skewed one- 
point-distribution function. We can also use the difference in the 
shape of the correlation function, since more extended correlation 
in the SZ effect differs from the point sources. The SZ power spec- 
trum scales roughly asT 2 while that of the point sources is con- 
stant, so the two can be distinguished. In practice noise and beam 
smoothing may make this analysis more complicated. 

Within a few years, several experiments will likely measure 
the SZ effect with an arcminute angular resolution and the effec- 
tive area of several square degrees. This has motivated us to revisit 
the predictions for the SZ angular power spectrum. Many analytic 
calculations exist in the literature (Cole & Kaiser 1988; Makino & 
Suto 1993; Atrio-Barandela & Mucket 1999; Komatsu & Kitayama 
1999; Cooray 2000; Molnar & Birkinshaw 2000; Holder & Carl- 
strom 2001; Zhang & Pen 2001), as well as many hydrodynamic- 
simulation results (Scaramella, Cen, & Ostriker 1993; Persi et al. 
1995; Refregier et al. 2000; Seljak, Burwell, & Pen 2001; da Silva 
et al. 2000; da Silva et al. 2001a; da Silva et al. 2001b; Springel, 
White, & Hernquist 2001; Refregier & Teyssier 2000; Zhang, Pen, 
& Wang 2002). The comparison between different simulations is 
difficult, as the parameters used differ, and the resulting power 
spectrum is extremely sensitive to the chosen parameters. Since 
massive clusters are at the tail of the dark-matter mass function, 
the r.m.s. mass fluctuations within a 8 h^ 1 Mpc sphere, as, is the 
most important parameter among the cosmological parameters. Ko- 
matsu & Kitayama (1999) have shown that the matter density of the 
universe is less important than ag. 

Refregier et al. (2000) and Refregier & Teyssier (2000) have 
attempted to compare analytic predictions based upon the so-called 
"halo approach" (section 2) with their simulated power spectrum, 
finding a good agreement. This agreement is, however, partially 
built into their model; in their approach, the parameters such as up- 
per and lower mass cut-off are chosen to agree with the simulations. 
Another free parameter used in the analytic modeling that cannot 
be obtained from these models is the mass-temperature normaliza- 
tion. Although a virial relation has often been used in this context, 
the overall amplitude is treated as a free parameter, and cannot be 
justified without a more detailed model. In this paper we develop a 
more refined analytic model and show that it reproduces the simu- 
lation results without any free parameters, at least at the level of the 
simulations disagreeing among themselves. 

The main ingredients into the halo model are the dark-matter- 
halo mass function, the gas-density and the gas-temperature pro- 
files. For the latter, we use an analytic gas and temperature profile 
model developed in our previous paper (Komatsu & Seljak 2001). 
We have shown that it reproduces well the simulation results in the 
outer parts of the cluster which are relevant for the SZ effect. This 
model fixes the normalization between mass and temperature by 
requiring that the gas and the dark-matter profiles agree with each 
other outside the gas core. Hydrodynamic simulations (e.g., Frenk 
et al. 1999) have shown that this assumption holds accurately. As a 
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result of this assumption and hydrostatic equilibrium the gas tem- 
perature declines in the outer parts of the cluster, which is also seen 
in the simulations (although it remains controversial in X-ray obser- 
vations due to the faint levels of X-ray emission in the outer parts of 
the cluster). This model therefore reproduces the main features of 
the hydrodynamic simulations and can be viewed as the gas analog 
of the often-used NFW profile (Navarro, Frenk, & White 1997) for 
dark matter, at least outside the gas core where additional physics 
may play a role. 

The structure of this paper is as follows. In § 2 we outline the 
"halo approach" to calculating the SZ angular power spectrum. In 
§ 3 we derive the universal gas-density, gas-temperature, and gas- 
pressure profiles following the prescription described by Komatsu 
& Seljak (2001). In § 4 we predict the spectrum and compare it 
to the hydrodynamic simulations. We then investigate the depen- 
dence of the spectrum on various cosmological parameters. We also 
study the mass and redshift distribution of the spectrum. In § 5 we 
compute the covariance matrix of Ci by taking into account non- 
Gaussianity of the SZ fluctuations and compare these to the hydro- 
dynamic simulations. In § 6 we show how well we can determine as 
with realistic SZ observations using a likelihood analysis. We put 
constraints on a Sl m — erg plane using the recent CBI and BIMA data 
of the small-scale CMB fluctuations at I ~ 2000 - 6000 (Mason 
et al. 2002; Dawson et al. 2002). In § 7 we quantify the theoretical 
uncertainty in our predictions. Finally, we summarize our results in 
§ 8. The fiducial cosmological model in this paper is f2 m = 0.37, 
Q A = 0.63, n b = 0.05, h = 0.7, to = -1.0, n = 1.0, and 
(78 = 1.0. 



2 HALO APPROACH TO THE SZ ANGULAR POWER 
SPECTRUM 

To compute the angular power spectrum of the SZ effect we use the 
halo formalism (Cole & Kaiser 1988; Makino & Suto 1993; Atrio- 
Barandela & Miicket 1999; Komatsu & Kitayama 1999; Cooray 
2000; Molnar & Birkinshaw 2000; Holder & Carlstrom 2001). For 
the angular scales of interest here, I > 300, the one-halo Pois- 
son term dominates Ci even after the subtraction of local massive 
haloes (Komatsu & Kitayama 1999). We thus neglect the halo-halo 
correlation term throughout the paper. 
The expression for Ci is given by 

2 /" Zmax , dV /" Mmax ,„,dn(M,z) , l2 _ 

C=£j o ^- z \ u _ dAI ' ^ ' <D 

where g v is the spectral function of the SZ effect (—2 in the 
Rayleigh-Jeans limit) (Sunyaev & Zeldovich 1980), V(£) the co- 
moving volume of the universe at z per steradian, dn{M, z)/dM 
the comoving dark-matter-halo mass function, and y\ (M, z) the 2D 
Fourier transform of the projected Compton y-parameter. The virial 
mass M used here and its mass function are described in more de- 
tail in section 2. 1 . 

For the upper integration boundary of redshift it suffices to 
take .Zmax = 10 (figure 5). For the mass integration boundaries, 
M min and M max , we find that M min = 5 x 10 11 h^ 1 M Q and 
M max = 5 x 10 15 h^ 1 Mq suffice to get the integral to converge 
on all angular scales, whereas M m i n = 5x 10 12 h^ 1 Mq is suffi- 
cient for I < 10 4 (figure 6). Since the integrals converge the results 
are not subject to the specific choice of integration boundaries. 

The 2D Fourier transform of the projected Compton y- 
parameter, yi = yi(M, z), is given by 



4%r s /*°° 2 , ,sin(lx/l s ) 
yi = -jr- J dxx y 3D (x) , (2) 

where y3o{x) is the 3D radial profile of the Compton y-parameter 
(equation 7 in section 3). Note that y3n(x) has a dimension of 
(length) _1 , while yi is dimensionless. In equation (2) we have 
used a; as a scaled, non-dimensional radius, 

x = r/r B , (3) 

where r s = r s (M, z) is a scale radius which characterizes the 3D 
radial profile and the corresponding angular wave number is 

l s = d A /r s . (4) 

Here d A = d A (z) is the proper angular-diameter distance. The 
projection of the 3D profile onto the sky in equation (2) was done 
within the Limber's approximation. In section 3, we derive yso{x) 
using the universal gas and temperature profile. 

2.1 Dark-matter halo mass function 

For the dark-matter halo mass function dn(M,z)/dM in equa- 
tion (1) we use the universal mass function (Jenkins et al. 2001), 
which was derived from iV-body simulations. We should be care- 
ful about the definition of "mass" when using the mass function 
derived from iV-body simulations (e.g., White 2001; Seljak 2001; 
Hu & Kravtsov 2002). The mass of haloes in the simulations is 
usually not the virial mass M as defined in our analytic formu- 
lation (equation 1). Instead, the simulations usually use either the 
friend-of-friend (FOF) mass or the spherical overdensity (SO) mass 
to define the haloes. We use here the SO mass, as the FOF mass is 
difficult to interpret in the context of analytic formulation. 

We define the SO mass, Ms, as the mass within the spheri- 
cal region whose density is S times the critical density of the uni- 
verse at z p c (z) (note that sometimes one uses instead the mass 
within the spherical region whose density is 5 times the mean mass 
density of the universe at z, p m («) = flm{z)p c {z); the two are 
related through S — O m (z)i5). Since TV-body simulations give 
dn(M s , z)/dM s , we need to convert it into dn(M,z)/dM. We 
do this by calculating 

dn(M,z) _ dM s dn(Ms,z) 
dM ~~ ~dM dMs ' 

For a given virial mass, M, we calculate Ms with equation (14) 
of Komatsu & Seljak (2001), and also evaluate dMs/dM numeri- 
cally. 

We use equation (B3) of Jenkins et al. (2001), which explicitly 
uses 8 — 180O m (2) (mass density of haloes is 180 times the mean 
mass density of the universe at z), 

dn(M s ,z) _ ^ pc(O) dlno-' 1 
dMs ~ M j Ms dMs 

x0.301cxp (- (mo- -1 +0.64| 3 ' 82 ) , (6) 

where a — a(Ms,z) is the linear, r.m.s. mass fluctuations at a 
given redshift z within the top-hat filter. To compute a we use the 
BBKS transfer function (Bardeen et al. 1986) with the baryonic 
correction made by Sugiyama (1995). The present-day critical den- 
sity of the universe is p c (0) = 2.775 x 10 11 h 2 Mq Mpc~ 3 . Note 
that we have defined the mass function as the comoving number 
density. 

While the functional formula for the mass function in equa- 
tion (6) was tested for just one cosmological model (rCDM model) 
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by Jenkins et al. (2001), it is essentially a non-smoothed version of 
their universal mass function valid for a broad range of models. Hu 
& Kravtsov (2002) have shown that the formula agrees with their 
simulated mass function as well which uses a different cosmolog- 
ical model. In addition, Evrard et al. (2002) have shown that their 
mass function for the SO mass haloes with S = 200 agrees with the 
Jenkins et al.'s formula, once taking into account the mass defini- 
tion properly (see also related discussion in Seljak 2001 and Hu & 
Kravtsov 2002). There is thus increasing evidence for the accuracy 
of Jenkins et al.'s mass function. 



3 GAS-PRESSURE PROFILE OF HALOES 
3.1 Introduction 

The 3D Compton y-parameter profile, y3o{x), is given by a ther- 
mal gas-pressure profile, P gas (a;), through 



y 3 u(x) = 



m e cr m a c 
1.04 x 10~ 4 Mpc -1 



/2_+2X\ 
K3 + 5XJ 

Pgzs(x) 



Pgns(x) 



50 eV cm -3 



(7) 



where P e (x) is an electron-pressure profile, X = 0.76 the primor- 
dial hydrogen abundance, or the Thomson cross section, m e the 
electron mass, and c the speed of light. 

We may write P gas (x) with a gas-density profile, p gas (x), and 
a gas-temperature profile, T gas (:r), as 



Pgas(x) 



3 + 5X p gas (a;) 



k B T gaa (x) 



55.0 eV cm 



• 0) 



10 14 M Q Mpc- 



feT gas (a;) 



8 kcV 



(8) 



where m p is the proton mass, and ks the Boltzmann constant. 

The aim of this section is to derive P gas (x) following Komatsu 
& Seljak (2001). Pressure profile was one of the uncertainties in 
the previous work on the SZ power spectrum: most of the previous 
work uses a spherical-isothermal (3 model as a gas-density profile, 
with P fixed at 2/3 for simplicity (Makino & Suto 1993; Atrio- 
Barandela & Miicket 1999; Komatsu & Kitayama 1999; Molnar & 
Birkinshaw 2000; Holder & Carlstrom 2001); Cooray (2000) used 
an isothermal gas-density profile as predicted by Makino, Sasaki, 
& Suto (1998), which uses hydrostatic equilibrium between gas 
pressure and dark-matter potential, and is close to the (3 profile. 
The p profile does not, however, reproduce the simulation results 
in the outer parts of the haloes. The outer slope of the /3 model 
with j3 = 2/3 scales as r~ 2 , being significantly shallower than the 
gas profile in simulations which, as we have argued, is similar to 
the dark-matter profile and asymptotically scales as r . Since the 
SZ effect is sensitive to the outer parts of haloes, this can cause a 
significant error in the predictions. 

There is a related uncertainty caused by the shallow (r~ 2 ) pro- 
file: one has to cut out the extension of gas in haloes at an arbitrary 
radius, otherwise one predicts too much power on large angular 
scales. We find that the previous predictions are unstable against 
this cut-off radius. This radius, which is a free parameter, has been 
assumed to be the virial radius; however, for isothermal gas as- 
sumed in the previous work this is ad hoc, as neither temperature 
nor density profile shows the abrupt cut-off at the virial radius in 
hydrodynamic simulations (Bryan & Norman 1998; Eke, Navarro, 
& Frenk 1998; Frenk et al. 1999). Instead, both the temperature 



and the density decrease smoothly across the virial radius, and it is 
only at 2 — 3 times the virial radius that there seems to be an abrupt 
decrease of the temperature to the IGM value, the so-called shock 
radius (Bryan & Norman 1998). 

Once the decline in profiles relative to the /3 model is included, 
the contribution from the outer parts of haloes to the SZ power 
spectrum converges well within the shock radius, and is thus inde- 
pendent of the outer cut-off of the profile. This is because requiring 
that the gas-density profile traces the dark-matter-density profile 



leads to p g 



in the outer parts of the haloes. This by itself 



would give a logarithmically divergent gas-mass profile; however, 
since the temperature is also declining with radius, the resulting 
pressure profile is convergent. Hydrostatic equilibrium and the gas 
density tracing the dark-matter density make the gas temperature 
decrease with radius. 



3.2 Universal dark-matter-density profile 

The universal NFW dark-matter-density profile (Navarro, Frenk, & 
White 1997) is 

Ps 



Pdm(x) = 



(9) 



x{l + x) 2 ' 

where x = r/r s , r s a scale radius, and p s a scale density. While the 
exact shape in the inner parts is still uncertain (Moore et al. 1998; 
Jing & Suto 2000; Klypin et al. 2001), it does not affect the SZ 
effect, given the scales of interest here. 

To specify the NFW profile, we need to specify a scale radius, 
r s , as a function of mass and redshift. It is customary to specify the 
concentration parameter c instead of r s , 

-0.2 



c(M,z) 



r (M, 



,(M,z) 



10 

1 + 2 



M 



M»(0) 



(10) 



where r v i r (M, z) is the virial radius of haloes, and the last expres- 
sion follows from Seljak (2000) with redshift evolution found by 
Bullock et al. (2001). The "non-linear mass" at z = 0, M*(0), is a 
solution to a(M) — 8 C , where a(M) is the present-day r.m.s. mass 
fluctuations within the top-hat filter, M the virial mass, and S c the 
threshold overdensity of spherical collapse at z — (Lacey & Cole 
1993; Nakamura & Suto 1997). We calculate the virial radius using 
the spherical collapse model, 

1/3 

r (M,z) 



M 



(4^/3)A c (z)p c (z) 



(11) 



where A c (z) is a spherical overdensity of the virialized halo within 
r v ; r at z, in units of the critical density of the universe at z (Lacey 
& Cole 1993; Nakamura & Suto 1997). 

The exact form of c(M, z), equation (10), is somewhat uncer- 
tain. There are various fitting formulae in the literature (Navarro, 
Frenk, & White 1997; Eke, Navarro, & Steinmetz 2001; Bullock 
et al. 2001; Seljak 2000; Cooray, Hu, & Miralda-Escude 2000). In 
section 4 we study the effect of different concentration parameters 
on our predictions, in addition to using our fiducial concentration 
parameter dependence in equation (10). We find that the difference 
in concentration parameters does not affect our predictions, largely 
because the SZ effect is not sensitive to the central parts of haloes, 
where the profile shape makes most difference. 

3.3 Universal gas-pressure profile 

To derive a gas-pressure profile P gas (x) we make three assump- 
tions: (1) hydrostatic equilibrium between gas pressure and dark- 
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matter potential due to the universal dark-matter density profile, 
(2) gas density tracing the universal dark-matter density in the outer 
parts of haloes, and (3) a constant polytropic equation of state for 
gas, P gas oc Pg as . The assumption (3) is obviously just an approx- 
imation and can be justified only to the extent that it suffices to 
explain the observations or simulations. There is evidence for this 
assumption to fail within inner 5% of the virial radius, where X-ray 
observations indicate that temperature increases with radius (Allen, 
Schmidt, & Fabian 2001). Outside the core region, nevertheless, the 
constant polytropic index seems to suffice, at least in comparison 
to the simulations. 

From these three assumptions one can obtain universal gas- 
density, gas-temperature, and gas-pressure profiles as 



pgas (a;) : 

Pgas (a?) 

and 



(0)j/ gas (x) 



Pgas(x) = P g as (0)?^ (a;)- 



(12) 
(13) 

(14) 



Solving the hydrostatic equilibrium equation with the ansatz above, 
we find that the non-dimensional gas-density profile, y sas (x), has 
an analytic solution, and is given by (Komatsu & Seljak 2001) 

1/(7-1) 



(X) 



1 



1 



ln(l + x) 



where the coefficient B is 



B = 3rj 



ln(l + c) 



1 + c 



(15) 



(16) 



In appendix, we give exact formulae for the polytropic index, 7, 
and the mass-temperature normalization factor at the centre, r/(0). 
Here we provide instead useful fitting formulae: 

7 = 1.137 + 8.94 x 10~ 2 ln(c/5) - 3.68 x 10~ 3 (c - 5), (17) 



and 



?7(0) = 2.235 + 0.202(c - 5) - 1.16 x 10~ 3 (c - 5) 2 



(18) 



These fitting formulae are valid for 1 < c < 25. Since 7 > 1 and 
the density decreases with radius, the temperature also decreases 
with radius. 

From 77(0), we can determine the central temperature, 

Tg a s(0), as 

4 



r g as(o) 



7/(0) 



GrripM 



3 + 5X 3r v 



8.80 keV T)(0) 



M/(10 15 ftT 1 M ( 



r vir /(l ft- 1 Mpc) 



(19) 



The central gas density, p ga s(0), is determined by equation (21) 
below. We then compute the central gas pressure, P gas (0), by sub- 
stituting p g as(0) and Tgas(0) into equation (8). This completely 
specifies the gas-density, the gas-temperature, and the gas-pressure 
profiles of haloes and their normalizations. 

We determine the amount of gas in haloes by requiring that the 
gas density in the outer parts of haloes is fib / fim times the dark- 
matter density. Unlike the j3 profile this is self-consistent, as our 
gas-density profile agrees with the dark-matter profile outside the 
gas core. The virial radius thus has no special importance, and the 
matching can be performed over a broad range of radii without af- 
fecting the results. We should reduce the gas fraction by the amount 
that has been converted to stars, which is of order 10—20% for mas- 
sive clusters of interest here (Fukugita, Hogan, & Peebles 1998). 



We choose not to apply this correction, so our predictions for the 
power spectrum are somewhat overestimated and the correspond- 
ing amplitude constraints underestimated. This however translates 
into only a few percent correction to the amplitude of as. 
We normalize the gas density at the virial radius, 



Pgas(c) 



Pgas(0)l/ g as(c) = ^-pdm(c) 

fi h M c 2 



ln(l + c) 



fi m 47rrJ ir (1 + c) 2 
We use this equation to obtain the central gas density, 

s(0) = 



1 + c 



(20) 



7.96 x 10 13 M Mpc~ 3 



fi b ft 2 



xc 



M/(10 15 ftT 1 Mq) 
[r vir /(l ft- 1 Mpc)] 3 



c 2 (l + c) 2 



ln(l + c) 



1+c 



(21) 



4 PREDICTIONS FOR THE SZ ANGULAR POWER 
SPECTRUM 

In this section, we predict the SZ angular power spectrum, Ci , by 
integrating equation (1), and compare the predictions with hydro- 
dynamic simulations. We then study the sensitivity of Ci to various 
cosmological parameters. We also study mass and redshift distribu- 
tion of Ci , finding that the power spectrum primarily probes mas- 
sive haloes in a moderate redshift universe. 



4.1 Comparison with hydrodynamic simulations 

We begin by comparing our predictions for Ci with six hydrody- 
namic simulations (four independent codes), listed in table 1. Fig- 
ure 1 compares our analytic predictions with the simulated power 
spectra. 

The predictions generally agree with the simulations within a 
factor of two. We find the best agreement for the nominally highest- 
resolution (512 3 ) simulation run by Zhang, Pen, & Wang (2002) 
using mesh-based MMH code (bottom-left panel of figure 1). The 
same code, but for different cosmological parameters and at lower 
(256 3 ) resolution, agrees well below I < 10000, but lacks power 
above that (Seljak, Burwell, & Pen 2001) (top-left panel). Note that 
the number of realizations for this simulation was small, which ex- 
plains lack of power at low / compared to even lower (128 3 ) res- 
olution used by Refregier et al. (2000) for the same code and the 
same cosmological parameters. 

For the 256 3 simulation we have created more 2-d maps, cov- 
ering the effective area of 1000 square degrees. We find that with 
more realizations C; at low I increases, becoming into better agree- 
ment with Refregier et al. (2000). This indicates that the sampling 
variance is very large on large scales, but is significantly reduced 
above / ~ 2000 — 3000. With more realizations, one typically 
increases power on large scales, which explains some of the dif- 
ferences among different simulations (they differ in the effective 
simulated area significantly; see table 1). Even so most of the sim- 
ulations are too small to include sufficient number of very massive 
halos and for these multiple maps from a single simulation do not 
reduce the sampling variance. This can explain why the simulations 
are still below the analytic model at low I. At least for MMH code, 
we find that the numerical resolution seems to be important at high 



6 E. Komatsu & U. Seljak 




Figure 1. Comparison between the predicted SZ angular power spectra 
(solid) and the simulations (dashed; see table 1 for the meaning of the labels 
shown in each panel, and for the cosmological parameters used). The dotted 
lines indicate r.m.s. errors of the simulated power spectra. The bottom-right 
panel scales all the power spectra by <r| ■ The agreement between 

the scaled power spectra indicates that this combination of cosmological 
parameters controls the amplitude of C; . 



I, and it causes some of the discrepancy between our model and the 
MMH simulations. 

The comparison with one of the two SPH simulations, da 
Silva et al. (2001b) using HYDRA code, shows a good agree- 
ment in shape, while the GADGET code (Springel, White, & 
Hernquist 2001) seems to predict more small-scale power above 
/ ~ 5000 compared to our predictions. We find that the shape of 
Ci of Springel, White, & Hernquist (2001) is similar to Refregier 
& Teyssier (2000), which uses a mesh-based hydrodynamic code, 
RAMSES, so the disagreement is unlikely to be due to simply SPH 
versus mesh-based codes. We note that the SPH simulation pre- 
sented in Bond et al. (2002) also agrees in shape with Springel, 
White, & Hernquist (2001). 

A possible explanation for the discrepancy is that some simu- 
lations such as MMH have not yet converged for very high /, while 
at low / sampling variance seems to be the most likely explana- 
tion. The discrepancy between our analytic model and GADGET 
and RAMSES could be caused by failure of our model at high /: 
clusters at higher redshift are less regular with more merging and 
substructure, which boosts C\ at high One way to account for this 
in analytic model would be to include subhalos within halos, but to 
properly calibrate this one would need to use hydrodynamic simu- 
lations. Here we argue that around / ~ 3000 the theoretical model 
and the simulations seem to agree rather well, certainly within a 
factor of two or better. Any remaining discrepancy at low I could 
be due to insufficient sampling in simulations. Although a factor 
of two may sound like a large uncertainty in the Ci prediction, it 
translates into a rather small uncertainty in erg, as discussed below. 



Figure 2. Dependence of the SZ angular power spectrum on erg. From top 
to bottom, the lines indicate cr s = 1.2, 1.1, 1.05, 1.0, 0.95, 0.9, and 0.8, as 
shown in the figure. 



4.2 Dependence on cosmological parameters 

In previous work it has been argued that the SZ power spectrum 
is very sensitive to as (Komatsu & Kitayama 1999; Seljak, Bur- 
well, & Pen 2001; Zhang & Pen 2001; Zhang, Pen, & Wang 2002). 
Figure 2 shows this dependence within the analytic model. For ex- 
ample, a 40% change in as causes one order of magnitude change 
in C'i, while a 10% change in <rg causes a factor of two change. We 
find that Ci oc a\ provides a good description of this scaling (see 
also the bottom-right panel of figure 1). The strong dependence of 
Ci on ag indicates that even if the numerical simulations are uncer- 
tain at the level of a factor of two, this translates into less than 10% 
systematic uncertainty in erg. 

This property of the SZ power spectrum should be compared 
with the current uncertainty in a s from the local cluster abundance 
studies, for which the systematic uncertainty is larger than 10% 
because of the uncertainty in the mass-temperature relation (e.g., 
Pierpaoli, Scott, & White 2001; Seljak 2001; Viana, Nichol, & Lid- 
die 2002). Weak lensing (e.g., Van Waerbeke et al. 2002) and di- 
rect mass-temperature calibration may soon provide a more reli- 
able method to determine agf^ 5 combination. On the other hand, 
we show below that Ci is not sensitive to fi m in the range of current 
interest. This property can break the degeneracy between the two 
parameters. 

Figure 3 shows dependence of Ci on £l m assuming a flat uni- 
verse, i.e., J1a = 1 — ^m- While Ci oc Q m for fi m > 0.4, we find 
that it is almost independent of fi m for 0.15 < fi m < 0.4. This 
is because the comoving volume of the universe increases rapidly 
with Qm decreasing in this region, canceling out the effect of Sl m 
on the dark-matter mass function. As this range of fi m is favored 
by current observations, Ci is practically independent of f2 m . 

Figure 4 shows dependence of C\ on equation of state w, Q^, 
h, and n. While varying Q b or h, we have fixed Q. h h at a fiducial 
value, 0.035, since we know that Ci scales as (flbh) 2 . The figure 
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Figure 3. Dependence of the SZ angular power spectrum on Q m . From top 
to bottom, the lines indicate Cl m = 0.97, 0.77, 0.57, 0.47, 0.37, 0.27, and 
0.17, as shown in the figure. While varying m , we have assumed a flat 
universe, i.e., Q/± = 1 — f2 m , and w = —1.0. 



shows any residual dependence of Ci on Qh or h. We find that Ci 
depends upon w, fl h , h, or n very weakly compared to as (compare 
figure 4 with figure 2). The only exception is if w > —2/3, which 
is however not favored by current observations (Bean & Melchiorri 
2002). 

We find that the dependence of overall amplitude of Ci on cos- 
mological parameters for 0.15 < O m < 0.4 is well approximated 

by 

l(l + l)Cl _ oon „ T ^2_7/»bfe\ 2 



330 fiK a s 



(0.15 < On < 0.4) ,(22) 



2tt ~~ r ° V 0.035/ 

at the peak around I ~ 4000. As we have shown in figure 3, C\ 
is rather insensitive to f2 m over this region. A fit for fi m > 0.4 is 
given by 



Cl OC (fib/l) 2 f2mO"g' 5nn 



-0.9<r 8 



(n m > o.4) . 



(23) 



4.3 Redshift and mass distribution 

Figure 5 shows the redshift distribution of Ci for a given 
dlnQ_ z^JdM^l\UM,z)\ 2 



(24) 



dim Jd^fdM^MM,*)?- 

We find that haloes at z ~ 1 determine C\ at I ~ 3000 (angular 
scales around 3')- Haloes at z ~ 2 dominate Ci at / = 10000. 
Even haloes at z ~ 3 have a non-negligible contribution to C; for 
I > 10000. 

Figure 6 shows the mass distribution of C\ for a given I, 

2 

(25) 



dlnC ; _ Mjdz^^lMM,z^ 



dlnM J d ^/ dM i^fe(M, Z )| 2 - 

We find that massive haloes between 10 14 h^ 1 Mq and 
10 15 h' 1 Mq dominate Ci at / = 3000 with a peak at 3 x 



Figure 4. Dependence of the SZ angular power spectrum on the equation of 
state of the dark-energy component, w, the baryon density f2b, the Hubble 
constant in units of 100 km s~ 1 Mpc -1 , h, and the primordial power- 
spectrum slope n. Values of each parameter are shown in each panel of the 
figure. While varying or h, we have Q^h = 0.035 fixed. The panels for 
Sib an d h are to show any residual dependence of C; on these two param- 
eters. This figure is to be compared with figure 2, showing how insensitive 
Ci is to w, f!b, h, or n compared to erg. 



1 = 500 

1=1000 

1=3000 

1=5000 

1=10000 




Figure 5. Redshift distribution of C\ . We plot d In C; /d In z, equation (24), 
for a given I. From left to right it is shown I = 500, 1000, 3000, 5000, 
10000, and 20000, as indicated in the figure. 
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Figure 6. Mass distribution of (7j . We plot d In C; /d In M, equation (25), 
for a given I. From right to left it is shown I = 500, 1000, 3000, 5000, 
10000, and 20000, as indicated in the figure. 



10 14 h- 1 M Q . For I = 1000 the peak is at 5 x 10 14 h' 1 Mq, 
while for / = 10000 at 10 14 h' 1 Mq. The large-angular-scale Q 
is thus sensitive to the presence of very massive and rare haloes, 
which explains why the sampling variance is so large on large 
scales. Much smaller angular scales (< 1' or I > 10000) receive 
significant contribution from smaller mass halos, where uncertain 
physics of gas cooling, heating, or star formation makes predictions 
unreliable; however, for the angular scales of 2' — 4' the dominant 
contribution is from massive (10 14 - 10 15 h- 1 M Q ) clusters for 
which these uncertainties are small and which are still so abundant 
that the sampling variance is small. Fortunately, this is precisely the 
range targeted by the next generation CMB experiments. 



5 COVARIANCE MATRIX OF THE SZ ANGULAR 
POWER SPECTRUM 

To determine the cosmological parameters from the power spec- 
trum we need to know its covariance matrix. If fluctuations are 
Gaussian, then we can calculate the variance from the spectrum 
itself. If fluctuations are non-Gaussian, then the covariance matrix 
is no longer diagonal. In this case we need to know not only Ci, 
but also specific configurations of the angular trispectrum, the har- 
monic transform of the four-point correlation function, to estimate 
the covariance matrix. 

The SZ fluctuations observed in the hydrodynamic simula- 
tions are very non-Gaussian (Seljak, Burwell, & Pen 2001; Zhang, 
Pen, & Wang 2002; White, Henquist, & Springel 2001). The scatter 
among the C\ estimates is large, especially on large scales, and the 
estimates are strongly correlated. We must thus calculate the full 
covariance matrix, taking into account the non-Gaussianity. Here 
we will use both the halo approach (Cooray 2001), and the hydro- 
dynamic simulations (Seljak, Burwell, & Pen 2001), to develop a 
reliable method for covariance matrix calculation. 



The covariance matrix of Ci, M u /, is given by (Cooray 2001) 

M W = ((C? ta -C«) (C? ta -C,,)) 



= f -1 

J sky 



2 (CO 2 
21 + 1 



5w + 



Tw 

47T 



(26) 



where the angular bracket denotes the ensemble average, / s k y is 
a fraction of the sky covered by an experiment (/sky = 1 for all 
sky), C° bs the observationally measured power spectrum, and Ci 
the theoretically predicted power spectrum (equation 1). Note that 
(C° bs ) = C'i. The second term in the r.h.s. represents the trispec- 
trum contribution. 



5.1 Halo approach to the angular trispectrum 

In this section we calculate M n i, equation (26), using the halo 
approach, and compare it with the hydrodynamic simulations. An 
angular-trispectrum configuration is represented by a quadrilateral, 
and is characterized by four sides and one diagonal (Hu 2001). 
Among the trispectrum configurations, those which constitute two 
lines whose lengths are I and respectively, and have zero diag- 
onal, determine the power- spectrum covariance matrix. We denote 
these as Tw (equation 26). 

Using the halo approach one obtains for the Poisson term 
which dominates the trispectrum (Cooray 2001), 



4 



d y pmax d n(M,z) 

dz— / dM — K — '—t 

dz JM min dM 

2 l„~ /if ~M2 



x\y l {M,z)\ 2 \y v {M,z)Y 



(27) 



where all the quantities have the same meaning as in equation (1). 
It follows from equation (26) that (47r/ s k y ) _1 Ti ; / determines 
the non-diagonal terms of Mw- Figure 7 shows Tw divided by 
C'iCii. To make it represent a ratio of the "non-Gaussian error" 
(the second term in M H i of equation 26) to the "Gaussian er- 
ror" (the first term) we have plotted T H i /(CiCv) multiplied by 
yj{2l + 1)(2Z' + 1)/8tt, which is this ratio when I = I'. We find 
that the ratio decreases with I as Z~ 0,95 , or equivalently, T a oc 
/- 195 (Ci) 2 ,forZ^ 10 3 . 

Figure 7 also shows how Mw scales with I for a given In 
the limit of T ;i <§C (C;) 2 , the plotted quantity in figure 7 reduces to 
the correlation coefficient of C;, rw, given by 

M w 



rw 



y/M u M v 



8tt 



CiCv 



(T u « Cf) . (28) 



As Tw /(CiCii) varies slowly with I for a given I', Ci and CV 
remain highly correlated even for a large value of \l — l'\ (Cooray 
2001). Because of this reason binning does not reduce the error bars 
as much as it would for Gaussian fluctuations. 



5.2 Comparison with the simulated error bars 

To compare the predicted error bars of Ci with the simulated ones 
one must take into account binning or band averaging that the simu- 
lations use to produce C\ . If the binned Ci at a given / is the average 
of Ci between I — Al/2 < I < l+Al/2, then the covariance matrix 
is 



MW » / S ky 



2(c ( ; 



(21 + l)Al 



Sw + 



Tw 
4tt 



(29) 
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Figure 7. The SZ an- 

gular trispectrum, T w , multiplied by {21 + 1){21' + l)/CiC t /. From 
top to bottom it is shown I' = 1000, 2000, 3000, 5000, 10000, and 20000, 
as indicated in the figure. The solid line shows oc i — °- 95 , which means that 
the ratio of the non-Gaussian error to the Gaussian error decreases with I as 
;~ ' 95 , and Tu scales as oc Z -1 - 95 (C() 2 . The plotted quantity reduces to 
the correlation coefficient of C; in the limit of T n -C (C;) 2 . 

Here, we have assumed that Al <C I. Since both C; and Tip are 
sufficiently smooth, the binning does not affect the values of C\ 
and T w as long as Al < I. We find that Al/l < 0.1 suffices for 
this approximation to hold. 

Figure 8 compares the predicted fractional errors of Ci, 
ACi/Ci = y/Tvhi/Ci, with the simulated ones. We find that the 
predictions agree with the simulations within a factor of two, sim- 
ilar to the agreement we have found for the spectrum itself. There 
is no clear trend in the residuals and we both ovepredict and under- 
predict the variance compared to the simulations. Since the vari- 
ance is very sensitive to the sampling errors, it is possible that some 
of the discrepancy is caused by insufficient sampling in the simu- 
lations. Figure 8 also plots Gaussian approximation to the errors, 
showing that our predictions give a significant improvement over 
the Gaussian approximation. We have also compared the predicted 
correlation coefficient, r^/ (equation 28), with the simulations of 
Seljak et al. (2001), finding a good agreement between the two and 
confirming that for high / the neighboring bins in C; are strongly 
correlated if the binning is sufficiently broad for the trispectrum 
term to dominate. Proper modeling of cross-correlations between 
the bins is essential for the parameter-error estimation and for the 
proper treatment of the statistical significance of results. 



6 DETERMINATION OF POWER-SPECTRUM 
AMPLITUDE 

In this section, we analyze how well we can determine as by mea- 
suring the SZ angular power spectrum with realistic observations, 
and apply it to the recent CBI results (Mason et al. 2002). 



Figure 8. Predicted fractional r.m.s. errors of Cj, ACi/Ci = y/Mu/Ci, 
in comparison with the simulated errors (see table 1 for the meaning of the 
labels shown in each panel, and for the cosmological parameters used). We 
have multiplied the simulated errors by the square-root of a fraction of the 
sky covered by the simulations, / B k y . The solid lines show the predicted 
fractional errors, while the dashed lines show the simulated ones. The filled 
circles show Gaussian approximation to the fractional errors, ~ (ZAi) -1 / 2 
for I 2> 1, where Al is the binning size shown in each panel. The predicted 
errors agree with the simulated ones within a factor of two, providing a 
significant improvement over the Gaussian approximation. 

We perform a least-square fitting to Ci by calculating x 2 , 

x ^Y,(C l -C l )(M-') ul (C l ,-C l ,), (30) 

i<v 

where Ci is a reference power spectrum which uses the fiducial 
cosmological parameters, or the actual CBI data. We determine the 
confidence levels of the parameter estimation by calculating Ax 2 - 

We vary r! m and erg, i.e., \ 2 = X 2 (^m, erg). We fix the other 
parameters at the fiducial values, Qa = 1 — f2 m , f2b = 0.05, 
h = 0.7, w — —1.0, and n = 1.0. For a given set of cosmological 
parameters, we calculate Ci and M u t and then obtain \ 2 - By ap- 
proximating probability distribution function of Ci with a Gaussian 
we have Ax 2 = 2.30, 6.17, and 11.8 for 68.3%, 95.4%, and 99.7% 
confidence levels (C.L.), respectively. 

First, we assume that experiments are sampling-variance lim- 
ited for 2000 < I < 3000, or 2000 < I < 5000; thus, we consider 
a few arcminutes angular resolution and ignore instrumental noise. 
To save computational time we do not evaluate equation (30) at 
every I, but at 10 points for 2000 < I < 3000, or 20 points for 
2000 < I < 5000. This choice gives Al = 111 and 158, respec- 
tively, in equation (29). Both Ci and T w are smooth so that this is 
a good approximation. 

For the sky coverage we consider 1 deg 2 or 10 deg 2 survey 
of the sky. These parameters are characteristic for the on-going or 
forthcoming experiments: CBI (Padin et al. 2001) for 2000 < I < 
3000 and 1 deg 2 , SZA (Holder & Carlstrom 2001) for 2000 < 
I < 3000 and 10 deg 2 , AMIBA (Zhang, Pen, & Wang 2002) or 
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Figure 9. Contours showing how well we can determine erg against S7 m by 
measuring the SZ angular power spectrum with realistic observations. We 
consider those observations which measure C; between 2000 < I < 3000, 
and survey 1 deg (top panel), or 10 deg (bottom panel) of the sky. The 
top panel is similar to the CBI experiment and the bottom panel to the SZA 
experiment. The solid line shows A\ 2 = 2.30 (68.3% C.L.), the dashed 
line Ax 2 = 6.17 (95.4% C.L.), and the dotted line Ax 2 = 11.8 (99.7% 
C.L.). The cross marks the fiducial cosmological model. We find that the 
contour is very sensitive to crs, yet almost independent of f2 m at around 
H m ~ 0.2 (see also section 4). We find that a CBI-type observation should 
be able to determine as to within 10% at 2<r, while a SZA-type to within 
5% at the same level. 

ACT (Lyman Page, private communication) for 2000 < / < 5000 
and 10 deg 2 . AMIBA and ACT experiments are actually going to 
survey about 100 deg 2 of the sky, so that the actual constraints to 
be obtained from these experiments would be better than we show 
here. 

Figure 9 shows \ 2 on a — fs plane for the 2000 < / < 
3000 region. The solid line shows Ax 2 = 2.30 (68.3% C.L.), 
the dashed line Ax 2 = 6.17 (95.4% C.L.), and the dotted line 
Ax 2 = 11.8 (99.7% C.L.). The top panel uses 1 deg 2 of the sky, 
close to the CBI experiment (if it is sampling-variance limited), 
while the bottom panel uses 10 deg 2 of the sky, close to the SZA 
experiment. We find that the contours are very sensitive to as and 
almost independent of £7 m for 0.15 < O m < 0.4 (Komatsu & 
Kitayama 1999). This is consistent with what we have previously 
shown in section 4 (see equation 22). While for one degree sky 
coverage the sampling-variance error is similar to the systematic 
(theoretical) error (10% at the 95% confidence level), a 10-times 
larger survey like SZA with the same angular resolution should be 
able to determine as to within 5% at the same level, so at that point 
the errors will be dominated by the theoretical uncertainty rather 
than by statistics. 

Figure 10 shows constraints on a Q m — as plane which would 
be obtained with ACT or AMIBA, whose angular resolution is 
1' - 2'. The figure uses 2000 < I < 5000 with 1 deg 2 (top panel) 
or 10 deg 2 (bottom panel) sky coverage. We find that the statistical 
error on as is extremely small for a 10 deg 2 survey. Surveys such 



Figure 10. The same figure as figure 9, but for 2000 < / < 5000. We find 
that statistical errors on as are very small for a 10 deg 2 survey, showing 
that ACT or AMIBA experiments which survey about 100 deg 2 of the sky 
determine as with theoretical-uncertainty-limited accuracy (~ 10%). A 
1 deg 2 survey shown in the top panel is sufficient enough to determine as 
with that accuracy. 



as ACT or AMIBA covering about 100 deg of the sky will deter- 
mine £78 with accuracy limited by the theoretical uncertainty. This 
is currently at a ~ 10% level and obviously more theoretical work 
should improve it in order for the power of the upcoming CMB ex- 
periments to be maximally exploited. It is of course still useful to 
cover such a large area of the sky, since this will allow us to use a 
bootstrap error determination without relying on the simulations. 

Recently, CBI and BIMA announced detections of CMB 
power at I ~ 2000 — 6000, which at face value are inconsistent 
with the primary CMB alone (Mason et al. 2002; Dawson et al. 
2002). Figure 11 compares the CBI and BIMA data (light-gray 
shaded area, which assumes Gaussian errors) to the predicted SZ 
angular power spectra for a s = 0.95, 1.05, and 1.15 (thin solid 
lines). We also show the primary CMB anisotropy for the fidu- 
cial cosmological model (dashed line), and the sum of the two for 
as = 1.05 (thick solid line). The predicted non-Gaussian errors are 
shown for as = 1.05 (dark-gray shaded area). We have included 
the instrumental-noise power spectrum in the non-Gaussian errors 
by using the noise power spectrum of 1(1 + l)C'i/(2n) = 500, 
1000, 1500, and 3000 fiK 2 for the CBI bins, and 720 and 2000 fiK 2 
for the BIMA bins. We have assumed the sky coverage of 1 deg 2 
for CBI and 0.1 deg 2 for BIMA. From the figure one can see that 
the data favor £78 ~ 1 models, and the predicted power spectra can 
explain the detected excess power reasonably well. 

We fitted the data to the predicted SZ angular power spec- 
trum using the four bins shown in figure 11 and computing x 2 
from equation (30). We compute the covariance matrix using equa- 
tion (29) with the bin-size of Al = 565, 378, 612, 1000, 2870, and 
4150 for each bin. For example, we have obtained r.m.s. error of 
327, 398, 373, 425, 370, and 299 fi K 2 at each bin for £7 8 = 1.05. 
The cross-correlation between the bins is most significant for the 
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Figure 11. The CBI (Mason et al. 2002) and BIMA (Dawson et al. 2002) 
data with Gaussian errors (light-gray shaded area), the predicted SZ angular 
power spectra for erg = 0.95, 1.05, and 1.15 (thin solid lines), the primary 
CMB anisotropy for the fiducial cosmological model (dashed line), and the 
sum of the two for erg = 1-05 (thick solid line). The dark-shaded area gives 
non-Gaussian errors predicted for erg = 1.05. 



first bin, for which the correlation coefficients with the other 3 CBI 
bins are 0.46, 0.41, and 0.29. The remaining 3 coefficients are be- 
low 0.3. Figure 12 shows confidence-level contours on the Sl m — og 
plane for the CBI and BIMA data. The top panel (a) assumes that 
the detection of the excess power in the all bins is entirely due to 
the SZ effect (i.e., we ignore the primary CMB anisotropy). We 
find CTg(fi b V0.035) 2/7 = 1-04 ± 0.12 at 95% confidence limit 
for 0.1 < a m < 0.5. 

The cluster abundance and cosmic-shear constraints are also 
shown in the figure as asQ.^ 5 — 0.45 ± 0.1, which roughly sum- 
marizes the current (systematic dominated) uncertainty from these 
methods. Combining the two can break the degeneracy between ug 
and fi m and is leading to cr 8 ~ 1, fi m ~ 0.2. One should be cau- 
tious not to overinterpret this conclusion, since systematic effects 
in both are still very large and could affect the parameter determi- 
nation. 

The bottom panel (b) assumes that the / ~ 1700 bin is af- 
fected by the primary CMB anisotropy, while the other bins are 
entirely due to the SZ effect. We assume that the primary CMB 
is 400 /iK 2 at the first bin, roughly consistent with the best-fitting 
theoretical curve, and is negligible in the other three bins. In prin- 
ciple we should change the primary CMB as well as we change 
the cosmological parameters; however, for simplicity we just use 
400 /iK 2 regardless of the cosmological parameters. We find a sim- 
ilar constraint on the fi m — cr 8 plane to the panel (a). In general, 
the CBI and BIMA data are consistent with detection of the SZ 
fluctuations, and imply crg(f2 b ft/0.035) 2/7 = 1.04 ± 0.12, but 
a somewhat larger survey is needed to confirm this interpretation. 
The good news is that with a survey just a few times the survey area 
of CBI one should be able to determine ug with an accuracy corn- 



Figure 12. The CBI and BIMA constraints on the Sl m — erg plane. The top 
panel (a) assumes that all the data shown in figure 11 are entirely due to 
the SZ effect. The dashed line shows A% 2 = 6.17 (95.4% C.L.) while the 
dotted line Ax 2 = 11.8 (99.7% C.L.). The cluster abundance and cosmic- 
shear constraints are shown in the figure as crgQg; 5 = 0.45 ± 0.1, which 
roughly summarizes the current uncertainty from these methods. The bot- 
tom panel (b) assumes that the first bin of the CBI data is affected by the 
primary CMB anisotropy. In both cases, we find crg(f2 b /i/0.035) 2 ' /7 = 
1.04 ± 0.12 at the 95.4% confidence level for 0.1 < J7 m < 0.5. 

parable to or better than any other current observations and with a 
different, and perhaps smaller, systematic uncertainty. 



7 UNCERTAINTIES IN THE PREDICTIONS 

The predictions are very sensitive to the halo mass function, so 
we investigate how the prediction changes if using different for- 
mulae in the literature. We compare four different predictions for 
Ci computed with four different formulae for the mass function: 
(IB3) Jenkins et al.'s tCDM mass function that uses the SO mass 
with S = 180f2 m (z) (Jenkins et al. 2001), (J9) Jenkins et al.'s uni- 
versal mass function that uses the FOF mass (equation 9 of Jenkins 
et al. 2001) and includes smoothing which may systematically in- 
crease the mass function in the exponential tail (Hu & Kravtsov 
2002), (HV) Evrard et al.'s Hubble-volume mass function which 
uses the SO mass with 8 = 200 (Evrard et al. 2002), and (PS) Press 
& Schechter's mass function which uses the virial mass (Press & 
Schechter 1974) (this one is of interest just for historical compar- 
ison). As Jenkins et al. (2001) note that the FOF mass they use 
corresponds to the SO mass with S = 180f2 m (z), we use the SO 
mass to calculate the mass function of (J9). 

The top-left panel (a) of figure 13 compares Ci computed with 
the four different mass functions. We find a very good agreement 
between (JB3) and (HV), while (J9) predicts higher amplitude. Hu 
& Kravtsov (2002) have also found that the mass function of (J9) 
predicts too many high-mass haloes to agree with their iV-body 
simulations, and concluded that the smoothing used for deriving 
(J9) causes the high-mass halo abundance to increase anomalously. 



12 



E. Komatsu & U. Seljak 



10" 10 



^3 

I 



\ 

T 1 



10" 



io- 




io- 



(a) mass function 



" — J(B3) PS 

--J(9) — -H.V. 

Hfmj — i i M i n i — m+ 




(c) gas extention 




(b) concentration 



S 



— B 



(d) excise inner 
region 



Figure 13. Theoretical uncertainty in the SZ angular power spectrum, (a) 
Uncertainty in the dark-matter-halo mass function. The solid line shows 
the prediction with Jenkins et al.'s tCDM mass function (equation 6), the 
dashed line Jenkins et al.'s universal mass function (equation 9 of Jenkins 
et al. (2001)), and the dot-dashed line Evrard et al.'s Hubble-volume mass 
function (Evrard et al. 2002). For comparison, the dotted line plots the pre- 
diction with Press & Schechter's mass function (Press & Schechter 1974). 
(b) Uncertainty in the concentration parameter, (c) Uncertainty in the hot- 
gas extension in haloes. The solid line shows C; with the hot-gas extension 
cut out at 3 times the virial radius, the dashed line at 2 times the virial ra- 
dius, and the dotted line at the virial radius, (d) Uncertainty in the inner gas- 
pressure profile. The solid line shows C; with no excision, the dashed line 
shows the one with the region interior to 10% of the virial radius excised, 
the dotted line 20% of the virial radius, and the dot-dashed line 40% of the 
virial radius. The dominant contribution to the spectrum is from 20 — 40% 
of the virial radius 



Jenkins et al. (2001) have explained this in their paper: when deriv- 
ing (J9), they are mainly interested in small-mass haloes for which 
the smoothing is less significant. Press & Schechter's mass func- 
tion predicts lower amplitude and underestimates the abundance 
of high-mass haloes compared with iV-body simulations (Jenkins 
et al. 2001). As the SZ effect is dominated by the high-mass haloes, 
we should not use (J9) and (PS), but (JB3) or (HV). 

Since there is no single expression for the concentration pa- 
rameter c (see equation 10 for definition) we need to quantify 
how much uncertainty in Ci is due to the uncertainty in the con- 
centration parameter. To quantify this, we use four different for- 
mulae for the concentration parameter: (S+B) power- spectrum- 
determined concentration parameter at z = (Seljak 2000) but 
evolved by (1 + z) _1 , for which evolution has been measured in 
TV-body simulations (Bullock et al. 2001), (S) power- spectrum de- 
termination with the non-linear mass evaluated at z, (B) halo de- 
termination of the concentration parameter (Bullock et al. 2001) 
and (CHM) a different power-spectrum-determined concentration 
parameter (Cooray, Hu, & Miralda-Escude 2000). The top-right 
panel (b) of figure 13 shows the predictions computed with the four 
different concentration parameters. We find that the uncertainty in 
the concentration parameter has little effect on Ci . 



How sensitive are the results to the outer radius of the gas pro- 
file (shock radius), where the gas temperature rapidly drops to IGM 
values? Previous work on a gas profile has assumed that the hot gas 
in haloes extends only up to the virial radius, but numerical sim- 
ulations suggest a smooth decrease in temperature at least out to 
twice the virial radius (Bryan & Norman 1998; Eke, Navarro, & 
Frenk 1998; Frenk et al. 1999). The bottom-left panel (c) of fig- 
ure 13 shows the predicted Ci with the hot-gas extension cut out at 
1, 2, or 3 times the virial radius. While a cut-off at the virial radius 
suppresses Ci at I < 5000, there is little difference for the cut-off 
at 2r v ir and beyond. One should thus include the hot-gas extension 
beyond the virial radius; otherwise, one underestimates the power 
spectrum on large angular scales, but the actual position of the cut- 
off is not very important since the power spectrum converges at 

2f vir • 

The gas density profile in the inner region displays a core 
(Bryan & Norman 1998; Eke, Navarro, & Frenk 1998; Frenk et al. 
1999), which is well reproduced by our model (Komatsu & Sel- 
jak 2001). On the other hand, recent X-ray observations have ob- 
served a decrease in temperature toward the centre within 5% of 
the virial radius (McNamara et al. 2000; Allen, Ettori, & Fabian 
2001; Allen et al. 2001; Arnaud et al. 2001a; Arnaud et al. 2001b; 
Kaastra et al. 2001), which our polytropic model does not account 
for. The bottom-right panel (d) of figure 13 shows the predicted C\ 
with the inner region within 10%, 20%, or 40% of the virial ra- 
dius excised. We find that the region interior to 10% of the virial 
radius has a negligible contribution to C\ , so the pressure profile in 
the cluster core has no effect on C; . From this figure we find that 
the dominant contribution to Ci comes from around 20 — 40% of 
the virial radius, roughly corresponding to r2soo, which is just in 
the range of Chandra satellite for medium redshift clusters (McNa- 
mara et al. 2000; Allen, Ettori, & Fabian 2001; Allen et al. 2001) 
and well within XMM/Newton or BeppoSax (Arnaud et al. 2001a; 
Arnaud et al. 2001b; Kaastra et al. 2001). Recent observations con- 
firm that the gas profile agrees well with the dark-matter profile in 
this region (Allen, Schmidt, & Fabian 2002), confirming our ba- 
sic assumption. We note that the temperature gradient is still small 
over this range, so most of the pressure gradient is caused by the 
gas-density profile. 

The effect of non-adiabatic physics has been investigated 
with hydrodynamic simulations. In da Silva et al. (2001b) it has 
been claimed that radiative cooling reduces Ci by modest amount, 
20 — 40%, on all angular scales. This is presumably due to a de- 
crease in the gas fraction, since some fraction of gas is being trans- 
formed into stars. Despite their extreme cooling treatment (giving 
too much cooled material to match the observations) the effect of 
radiative cooling is modest (note that 40% change in Ci gives a 5% 
change in erg). Preheating may lead to a larger effect, as it may al- 
ter the gas structure not only in the inner region, but also at larger 
radii, especially for smaller-mass haloes. Some fraction of the gas 
may be expelled from the cluster, which would suppress the SZ 
fluctuations on small scales. At the same time, preheating increases 
the gas temperature, which would increase the spectrum. However, 
care must be taken not to violate the FIRAS constraints on the mean 
y parameter and it is unlikely that this could significantly enhance 
the SZ power spectrum without violating these limits. Several au- 
thors have noticed significant effects of preheating on Ci (Springel, 
White, & Hernquist 2001; Holder & Carlstrom 2001; da Silva et al. 
2001b), although current treatments are still over-simplified. One 
must be careful in analytic models to make them self-consistent, 
as one cannot retain gas in hydrostatic equilibrium within exter- 
nal dark-matter potential using a large gas core and isothermal gas. 
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When hydrostatic equilibrium is imposed the pressure profile is less 
affected than the density or temperature profiles. It is likely that 
any preheating will have only a modest effect on the SZ power 
spectrum. Recently, White, Henquist, & Springel (2001) have used 
hydrodynamic simulations to show that gas cooling, energy feed- 
back, and star formation affect Ci by no more than a factor of two, 
confirming our conclusion. 



8 CONCLUSIONS 

In this paper we present a refined analytic model for the angular 
power spectrum of the SZ effect. In contrast to previous model- 
ing (Atrio-Barandela & Miicket 1999; Komatsu & Kitayama 1999; 
Cooray 2000; Molnar & Birkinshaw 2000; Holder & Carlstrom 
2001; Zhang & Pen 2001), our model has no free parameters and 
treats the gas in a cluster dark-matter potential self-consistently, us- 
ing the universal gas-density and temperature profile (Komatsu & 
Seljak 2001), which fits the observed and the simulated properties 
of clusters of galaxies in the range of interest to the SZ observations 
(outside 5% of the virial radius). 

We compare our predictions for the spectrum and the errors 
with the hydrodynamic SZ simulations and find a good agreement. 
The deviations between our model and the simulations are com- 
parable to the deviations among the simulations themselves and 
are particularly small in the observationally relevant range around 
/ ~ 3000. Some of the discrepancies are due to either poor reso- 
lution in the simulations or sampling variance, while the remain- 
ing difference could be due to the simplifying assumptions of gas 
physics in our model, such as hydrostatic equilibrium, spherical 
symmetry, merging and substructure, which may be inaccurate in 
particular at higher redshifts when the clusters are still forming. 
The disagreement is more important on very small angular scales 
(I > 5000), while the accuracy of our predictions in the range of 
I ~ 2000 — 5000 is better than a factor of two, which translates 
to less than 10% systematic uncertainty in the amplitude of fluc- 
tuations as- We also compare analytic predictions for the power- 
spectrum covariance matrix to the simulations and find similarly 
good agreement. 

We investigate the dependence of Ci on various cosmological 
parameters, finding that over the range of interest as and Q b h de- 
termine the amplitude of Ci almost entirely, as has been pointed 
out by several authors before as well (Komatsu & Kitayama 1999; 
Seljak et al. 2001; Zhang & Pen 2001; Zhang et al. 2002). It is par- 
ticularly important that the matter density of the universe does not 
affect the spectrum as much as as does, in agreement with Komatsu 
& Kitayama (1999). This differs from the local-cluster-abundance 
studies, where the constraint is usually on crgf^ 5 . 

The dominant contribution to Ci comes from massive clus- 
ters (M > 10 14 h' 1 M Q for / < 5000) at moderate (z ~ 1) 
redshift. Within the cluster, the dominant contribution comes from 
around 0.2 — 0.4r v i r , suggesting that the physical processes tak- 
ing place in the cluster core such as cooling, heating or heat con- 
duction have little effect on the SZ power spectrum. We assume 
that the gas fraction is close to the cosmic mean baryon fraction, 
which is reasonable for massive clusters and has observational sup- 
port in low stellar-mass-to-cluster-mass ratios in massive clusters 
(e.g., Fukugita, Hogan, & Peebles 1998), as well as in studies of 
the gas fraction as a function of radius (e.g., Allen, Schmidt, & 
Fabian 2002). While more exotic mechanisms such as preheating 
could affect C; on all scales (by increasing the temperature and/or 
reducing the gas clumping), they are constrained by the FIRAS lim- 



its and are unlikely to make more than a factor of two change in Ci 
(comparable to a 10% change in as). 

Using analytic predictions for the spectrum and its covariance 
matrix we have performed a likelihood analysis to estimate how 
well we can measure as against Sl m with realistic SZ observations. 
The likelihood is very sensitive to as almost independent of Q m 
(figure 9 or 10). By performing the likelihood analysis on the re- 
cently reported CBI (Mason et al. 2002) and BIMA (Dawson et al. 
2002) detections at the level of 15 - 20 fiK at I ~ 2000 - 6000 
(figure 11), we find as(fl b h/0.035) 2/7 = 1.04 ± 0.12 at the 95% 
statistical confidence level. To this we should add about 10% sys- 
tematic theoretical uncertainty due to the simplified modeling in 
our model, the numerical issues in the simulations, and missing 
physics in both. We should note that most of the effects we have ig- 
nored within our model tend to further increase as'- assumed fidu- 
cial value for £l b h = 0.035 is somewhat high, we ignored fraction 
of gas transformed into stars, and our predictions are somewhat 
higher than the simulation results at the angular scales relevant for 
CBI and BIMA. These effects would further increase as, although 
not by more than 5 — 10%. Such a high value for as is in tension 
with the primary CMB amplitude determination (Lahav et al. 2002) 
and Ly-a forest determination (Croft et al. 1999; McDonald et al. 
2000), although it can be accommodated by the cluster abundance 
and the weak lensing results if fi m ~ 0.2 ± 0.1. While the CBI 
experiment is surveying 1 deg 2 of the sky and still has a large in- 
strumental noise, a somewhat larger sky coverage with lower noise 
(e.g., SZA) would significantly reduce the statistical errors on as- 
Surveys with 100 deg 2 such as ACT or AMIBA should be able to 
measure both the power spectrum and its covariance matrix (e.g. by 
using bootstrap sampling) to exquisite precision. 

In conclusion, measurement of the angular power spectrum of 
the SZ effect offers a promising way to determine as, since it is 
very sensitive to its value and is almost independent of fi m . It is 
free of observational selection effects such as flux, surface bright- 
ness, or volume limit of the survey. It does not require the mass 
of sampled haloes to be measured, and there is no uncertainty as- 
sociated with the mass-flux calibration. Finally, it is within reach 
of the next generation of small-scale CMB experiments and may 
have already been detected by the CBI and BIMA experiments, in 
which case their results imply that the amplitude of fluctuations at 
8 hT 1 Mpc scale is of order unity. 
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Table 1. Simulation parameters. The first column assigns a label to each simulation. The second column shows the name of a code which each simulation uses: 
MMH (=Moving-Mesh Hydrodynamic code) (Pen 1998), and RAMSES (adaptive mesh-refinement hydrodynamic code) (Teyssier 2002) are mesh codes, 
while GADGET (=GAlaxies with Dark matter and Gas intEracT code) (Springel, Yoshida, & White 2001), and HYDRA (Pearce & Couchman 1997) are 
SPH (=Smoothed-Particle Hydrodynamic) codes. The third column shows the number of grids for mesh codes (MMH and RAMSES), or the number of 
SPH particles for SPH codes (GADGET and HYDRA). The box size of simulations is 100 h^ 1 Mpc except for SWH01, which uses 134 h^ 1 Mpc. The 
fourth column shows the number of two-dimensional maps which each simulation has created. The fifth column shows the field-of-view of the simulated 
two-dimensional maps. Note that RKSP00 and RT00 have not created two-dimensional maps, but computed the SZ angular power spectrum from the three- 
dimensional SZ power spectrum. The sixth— eighth columns show the cosmological parameters which each simulation uses. All the simulations use fl\ = 
1 — f2 m , and n = 1.0. The rightmost column shows references for the simulations. 



Label 


Code 


/Ygrid or Af sph 


■^V map 


F.O.V. 


f2 m 


n b 


h 


CT8 


Reference 


RKSP00 


MMH 


128 3 






0.37 


0.049 


0.7 


0.8 


Refregier et al. (2000) 


SBP01 


MMH 


256 3 


12 


2° X 2° 


0.37 


0.049 


0.7 


0.8 


Seljak, Burwell, & Pen (2001) 


SWH01 


GADGET 


224 3 


15 


1° x 1° 


0.3 


0.04 


0.67 


0.9 


Springel, White, & Hernquist (2001) 


daSilvaOl 


HYDRA 


160 3 


30 


1° x 1° 


0.35 


0.0377 


0.71 


0.9 


da Silva et al. (2001b) 


RT00 


RAMSES 


256 3 






0.3 


0.039 


0.7 


0.93 


Refregier & Teyssier (2000) 


ZPW02 


MMH 


512 3 


40 


1°19 x 1°19 


0.37 


0.05 


0.7 


1.0 


Zhang, Pen, & Wang (2002) 
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APPENDIX A: GAS POLYTROPIC INDEX 



In this appendix we derive an exact formula for the gas polytropic 
index 7 and give fitting formulae for 7 and for the mass-central 
temperature normalization factor, r)(0) (equation 19). 

We specify 7 and 77(C)) uniquely by requiring that the gas- 
density profile, p gas (x), matches the dark-matter density profile, 
pdm(x), in outer parts of the haloes. We do this by solving the fol- 
lowing equation, 



dlnp gas (a;) 



din a; 



din pdm(x) 



din a; 



(Al) 



where s» denotes a slope of the dark-matter density profile at a 
matching radius x*. 

In previous paper (Komatsu & Seljak 2001) we have shown 
that equation (Al) determines 77(0) uniquely for a given 7, 
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where m(x) is a non-dimensional dark-matter mass profile defined 
by 



m(x) 



duU Pdm(w), 



(A3) 



and p B is a scale density. Equation (A2) makes the matching condi- 
tion, equation (Al), satisfied at radius x*. 

Next, we require that 77(0) does not depend upon a particular 
choice of x*; this requirement specifies 7 uniquely. While we have 
done this by finding an empirical fitting formula for 7 in the previ- 
ous paper (equation 25 of Komatsu & Seljak 2001), we derive here 
an exact formula for 7 by solving the following equation, 
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= 0, 



(A4) 



where c is the concentration parameter (see equation 10 for defini- 
tion). This equation means that 77(0) normalization does not depend 
upon the particular choice of x*. 

We find that equation (A4) gives 7 as a function of c, 



din [m(x*)/s* 



s,dlnx* 



(A5) 



This formula gives the polytropic index of gas in haloes with any 
given dark-matter density profiles. 

Hereafter, we evaluate 7 and 77(0) for the NFW profile 
(Navarro, Frenk, & White 1997), p dm (x) = psX' 1 ^ + x)~ 2 . We 
find 

l + 3:r* 



and 
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(A6) 



(A7) 



By substituting s« and m(x) into equation (A5), we obtain 7. We 
then substitute s r , m(x), and the derived 7 into (A2) to obtain 77(0). 

The upper panel of figure Al shows 7 as a function of c 
compared to the fitting formula in equation (17). For comparison, 
the figure also plots the fitting formula which we have given for 
4 < c < 11 in Komatsu & Seljak (2001). The new formula, equa- 
tion (17), is valid for 1 < c < 25. 
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Figure Al. The top panel shows the polytropic index of gas in haloes, 7, 
while the bottom panel shows the mass-central temperature normalization 
factor, 77(0) (equation 19), as a function of the concentration parameter c. 
The solid lines plot exact values calculated from equation (A5) and (A2) 
for 7 and 77(0), respectively. The dotted lines plot the best-fitting formulae, 
equation (17) and (18) for 7 and 7/(0), respectively. The dashed lines plot 
the fitting formulae given in Komatsu & Seljak (2001) for 4 < c < 11. 



By substituting 7 from equation (A5) into equation (A2) we 
obtain 77(0) as a function of c. The bottom panel of figure Al shows 
77(0) as a function of c and also the fitting formula 77(0) in equation 
(18). For comparison, the figure also plots the fitting formula given 
in Komatsu & Seljak (2001). 



